fnames_in = {'hide\2018-02-01_11_56\expt_hydro_fit_all_four_params_bootstrap_summary_natural_units.txt',...
    'hide\2018-02-01_11_49\expt_hydro_fit_Gamma_k_sqrd_term_bootstrap_summary_natural_units.txt',...
    'hide\2018-02-01_11_46\expt_hydro_fit_2paramsonly_bootstrap_summary_natural_units.txt'};
%     'expt_hydro_fit_gammak_minlambda=0.0_summary_natural_units.txt'};

leg = {'\Gamma, D, for \lambda > 0',...
    '\Gamma, D, D_m fpr \lambda > 0'};
%     '\Gamma_o, \Gamma_1, D \lambda > 0'};

% create subplots and handles
figh = figure;
ax1 = subplot(2, 2, 1);
ax2 = subplot(2, 2, 2);
ax3 = subplot(2, 2, 3);
ax4 = subplot(2, 2, 4);


for ii = 1:length(fnames_in)
    fname_in = fnames_in{ii};
    dat = dlmread(fname_in, '\t', 1, 0);
    T = dat(:, 1);
    T_unc = dat(:, 2);
    gamma = dat(:, 5);
    gamma_unc = dat(:,6);
    d = dat(:, 7);
    d_unc = dat(:, 8);
    gamma_2 = dat(:, 9);
    gamma_2_unc = dat(:, 10);
    omega_2 = dat(:, 11);
    omega_2_unc = dat(:, 12);

    % plot data
    errorbar(ax1, T, gamma, gamma_unc, gamma_unc, T_unc, T_unc, 'o');
    hold(ax1, 'on');
    grid(ax1, 'on');
    ax1.YLim(1) = 0;
    ax1.YLabel.String = 'Gamma (t/\hbar)';
    ax1.XLabel.String = 'T (t)';

    errorbar(ax2, T, d, d_unc, d_unc, T_unc, T_unc, 'o');
    hold(ax2, 'on');
    grid(ax2, 'on');
    ax2.YLim(1) = 0;
    ax2.YLabel.String = 'D_ (ta^2/\hbar)';
    ax2.XLabel.String = 'T (t)';

    errorbar(ax3, T, gamma_2, gamma_2_unc, gamma_2_unc, T_unc, T_unc, 'o');
    hold(ax3, 'on');
    grid(ax3, 'on');
    ax3.YLim(1) = 0;
    ax3.YLabel.String = 'Gamma_2 (ta^2/\hbar)';
    ax3.XLabel.String = 'T (t)';

    errorbar(ax4, T, omega_2, omega_2_unc, omega_2_unc, T_unc, T_unc, 'o');
    hold(ax4, 'on');
    grid(ax4, 'on');
    ax4.YLim(1) = 0;
    ax4.YLabel.String = 'Omega_2 (ta^2/\hbar)';
    ax4.XLabel.String = 'T (t)';
end

legend(leg);


